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Abstract 

We examine possible J PC quantum number assignments for the X(3872). Angular correlations 
between final state particles in X(3872) — ► ir + ir~ J/ip decays are used to rule out J PC values of ++ 
and CT + . The shape of the tt + tt~ mass distribution near its upper kinematic limit favors 5-wave 
over P-wave as the relative orbital angular momentum between the final-state dipion and J/ip, 
which strongly disfavors 1 ^ and 2 h assignments. The accumulated evidence strongly favors a 
J PC = 1++ assignment for the X(3872), although the 2 ++ possibility is not ruled out by tests 
reported here. The analysis is based on a sample of X(3872) mesons produced via the exclusive 
process B — > KX(3872) in a 256 fb _1 data sample collected at the T(45) resonance in the Belle 
detector at the KEKB collider. 

PACS numbers: 14.40. Gx, 12.39.Mk, 13.20.He 
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The X(3872) was first observed by Belle in exclusive B~ —> K~7i + ir~ J/ip decays pj, 
The subsequent observation of the X(3872) — > jJ/ip decay mode established the 
charge parity as C = +1. In the same paper, Belle also reported evidence for the decay 
X — > 7t + 7t~tt°J/iI), where the 7r + 7r~7r° invariant mass distribution has a strong peak between 
750 MeV and the kinematic limit of 775 MeV, suggesting that the process is dominated by 
the sub-threshold decay X — > uJ/ip. The partial widths for "inJ/ip and 2irJ/%p decays are 
of comparable size, which implies a large violation of isospin symmetry. 

Here we report on a study of X(3872) — ► -n + n~J/ip decays produced via the exclusive 
decay process B —>■ KX (3872). We use a data sample that contains 275 million BB pairs 
collected in the Belle detector at the KEKB energy-asymmetric e + e~ collider. The data were 
accumulated at a center-of-mass system (cms) energy of \/s = 10.58 GeV, corresponding to 
the mass of the T(45) resonance. KEKB is described in detail in ref. |4j. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector, a 50-layer cylindrical drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect Kl mesons and to identify 
muons (KLM). The detector is described in detail elsewhere 

We select events that contain a J/tp-, either a charged or neutral kaon, and a 7r + 7r _ pair 
using criteria described in refs. [l[ and p. To reduce the level of e + e~ — > qq (q = u, d, s or c- 
quark) continuum events in the sample, we also require R2 < 0.4, where R2 is the normalized 
Fox- Wolfram moment 0, and | cos##| < 0.8, where Ob is the polar angle of the 5-meson 
direction in the cms. 

Candidate B — > Ktt + 7t~ J ' /ip mesons are identified by the energy difference AE = E B ms — 
-^beam an d the beam-energy constrained mass M bc = J (i^cam) 2 ~ (p^) 2 ! where E^ m is the 
cms beam energy, and E c ™ and p c ™ are the cms energy and momentum of the Kir + ix~ J /ip 
combination. We select events with M^ c > 5.20 GeV and \AE\ < 0.2 GeV and among 
these define a signal region 5.2725 GeV < M bc < 5.2875 GeV and \AE\ < 0.034 GeV; this 
corresponds to ±3<j from the central values for each variable. 

We select events with a dipion invariant mass requirement of M n + n - > (M (n + n~ J / ip) — 
(ttim + 200 MeV), which corresponds to M n + n - > 575 MeV for the A(3872). This reduces 
misidentified 7 conversions and combinatoric backgrounds by 36% with an V(3872) signal 
loss of 6%. 

These selection criteria isolate a very pure sample of 696 ±26 B — >• Ki/j(2S), ip{2S) — ► 
7T + 7T~J/if) events. These events are used as a calibration reaction to determine the Mb c , AE 
and M(n + 7T~ J/'ip) peak positions and resolution values, and for validating the Monte-Carlo 
(MC) acceptance calculations. 

Figure^shows the M(7r + 7r~ J / ip) mass distribution near 3872 MeV for the selected events. 
Here the smooth curve is the result of a fit with a Gaussian function to represent the 
V(3872) signal and a first-order polynomial to represent the background. The width of 
the Gaussian is fixed at a = 3.2 MeV, the experimental resolution determined from the 
ip(2S) — > 7T + 7T~J/if) event sample. The total signal yield is 49.1 ±8.4 events. For subsequent 
analysis, we define an A(3872) signal region to be ±5 MeV around the signal peak. For 
background estimates, we use ±50 MeV sidebands above and below the signal peak centered 
at 3837 MeV and 3907 MeV. There are a total 58 events in the signal region; the background 
content, determined from the scaled sidebands, is 11.4 ± 1.1 events. 
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FIG. 1: The M(vr + 7r" J/tp) mass distribution for the X(3872) region. 

Using a MC-determined acceptance, we determine the product branching fraction 
B(B -> XX(3872)) x B(X -> Tr+Tr" J/VO = 

1.31 ± 0.24(stat) ± 0.13(syst) x 1(T 5 . (1) 

where we have assumed equal B — > KX branching fractions for charged and neutral B 
mesons, and that the dipion originates from p — > 7r + 7r~. The systematic error includes the 
effect of uncertainties in the M(7r + 7r~) shape for X(3872) decay. This result agrees with, 
and supersedes, the results of ref. [l|. 

Since both the B and K mesons are scalar particles, X(3872) mesons produced via 
exclusive B — * KX decays cannot have a non-zero component of angular momentum along 
their momentum direction in the B rest frame. This provides useful limits on the number 
of independent partial- wave amplitudes needed to describe the decay pi. IqL E0|. 

With less than fifty signal events, any angular distribution will have, on average, only 
about five signal events per bin, which is not sufficient for a standard angular analysis. 
However, because the signal-to- noise ratio for the X — > it + -n~J/ip signal is quite good 
(S/N ~ 4), a typical distribution has, on average, only about one or two background events 
per bin. We exploit this good S/N and try to find, for a given J PC hypothesis for the 
X(3872), angular quantities that have distributions with a zero in some location. In the 
bins near the zero point, any observed events would have to be accounted for by upward 
fluctuations of the background [11 . 

For h , there is only one invariant amplitude corresponding to a p and J/if> in a P- 
wave. The decay amplitude is proportional to the scalar triple product of the p and J /if) 
polarizations and their relative momentum. As a result, the polarizations are perpendicular 
to each other and their relative momentum. We follow a suggestion by Rosner [^] and use a 
coordinate system where the x-axis is defined to be opposite the J /if) direction in the p rest 
frame, the x — y plane is defined by the ir + and J/if> directions and the z-axis is chosen so 
that it forms a right-handed coordinate system. We define 8 as the angle between the £ + 
and the z axis in the J/ if> rest frame and if> as the angle between the 7r + and the x axis in the 
dipion rest frame. The expected distribution for h is d 2 N/d(cos9)d(cosip) oc sin 2 9 sin 2 if). 

The | cos#| and | cos^l distributions for the X(3872) signal region are shown in Figs.|2fa) 
and (b), respectively. The shaded histograms indicate the side-band determined background. 
The distributions for both variables show strong signals at the upper edge of each plot, in 
contrast to expectations for a sin 2 6 sin 2 if) dependence. The open histogram shows the 0~ + 
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FIG. 2: The (a) |cos#| and (b) |cos?/>| distributions for events in the X(3872) signal region 
(points with error bars). The open histogram is the expected distribution for a h assignment 
including background. The hatched histogram shows the scaled sideband. 

MC expectations plus background, normalized to the observed number of events. Here the 
agreement is marginal for cos#: y^/d.o.f. = 17.7/9 but poor for cosip: x 2 /d-o.f. = 34.2/9. 
This latter distribution allows us to reject the 0~ + assignment with high confidence. 




0.00 0.50 1.00 

Icosej 

FIG. 3: The | cos distribution for events in the X(3872) signal region (points with error bars). 
The open histogram is the expected distribution for a ++ assignment including background. The 
hatched histogram shows the scaled sideband. 

For ++ , two invariant amplitudes are possible, corresponding to the p and J/ip in relative 
S- or .D-waves. Because of the limited phase-space, the D-wave contribution can be expected 
to be strongly suppressed relative to the S wave term and is ignored. The amplitude is then 
proportional to the scalar product of the p and J/ ip polarizations. We define 9e n as the angle 
between the £ + and the 7r + in the X(3872) rest frame. In the limit where the X(3872), J/ip 
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and p rest frames coincide dN/d(cos9i n ) oc sin 2 ^. The kinematic smearing due to relative 
motion of the different frames is incorporated in the MC simulations that are used to compare 
data with expectations |l3j |. 

Figure El shows the | cos^l distribution, computed in the p rest frame, for X(3872) signal 
region events. The agreement with S-wave ++ MC expectations is poor: x 2 /d.o.f. = 31.0/9, 
and provides evidence against the ++ assignment. 

For 1 ++ the J/ip and p can be in a relative S and/or D-w&ve. We use a coordinate 
system where the x-axis is the negative of the kaon flight path, the x — y plane is defined 
by the kaon and 7r + and the z axis completes a right-handed coordinate system. The angle 
between the 7r + direction and the x-axis is x an d the angle between the £ + direction and the 
z-axis is 9 1. In the limit where the J/ip and p are at rest in the X rest frame (and .D-wave con- 
tributions can be neglected), the amplitude is proportional to the vector triple product of the 
X, p and J/ip polarizations, and the choice of axes ensures that the X polarization is along 
the x direction 0, The expectation for 1 ++ is d 2 N / d{cos6^)d{cosx) oc sin 2 #^ sin 2 x. 




FIG. 4: The a) | cos^| and b) | cosxl distribution for events in the X(3872) signal region (points 
with error bars). The open histogram is the expected distribution for a 1 ++ assignment including 
background. The hatched histogram shows the scaled sideband. 

The | cos 9(\ distribution for X(3872) signal region events is shown in Fig.|3Ja). The distri- 
bution tends toward zero at the upper edge of the plot, as expected for a sin 2 $e dependence. 
The open histogram shows the results of a comparison to normalized MC expectations for 
1 ++ decaying to a p and J/ip in an S'-wave. The agreement is good: x 2 /d.o.f = 11.4/9. The 
| cosxl distribution is shown in Figs. IHb) together with the MC expectation for 1 ++ . The 
agreement here is also good: x 2 /d.o.f . = 5.0/9. 

For even-parity C — +1 states the 7r + 7r~ J/ip final state would be a p and J/ip primarily 
in a relative S'-wave, with some possible D-wave component. For odd-parity states the p 
and J/ijj would be in a relative P-wave with some possible F-wave. The M(7r + 7r~) mass 
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distribution near the upper kinematic boundary is suppressed by a (qj/^) 2£+1 centrifugal 
barrier, where (fji^ is the J/ip momentum in the X(3872) rest frame, and £ is the orbital 
angular momentum. For the S'-wave (i.e. J p = J + ) cases, the upper-boundary is modulated 
by the available phase-space, which is proportional to q*j/^] for a P-wave the modulation 
is (<2j/^) 3 - Thus, the shape of the high-mass part of the -n + n~ invariant mass distribution 
provides some J PC information. 




FIG. 5: M(ir + ir~) distribution for events in the X(3872) signal region; the histogram indicates 
the side-band determined background. The solid (dashed) curve shows the fit that uses a p Breit- 
Wigner line shape with the J ftp and p in a relative S- wave (P-wave). The dot-dashed curve is a 
smooth parameterization of the background that is used in the fit. 



Figure El shows the distribution for events in the X(3872) — > 7T + 7T~J/ifj signal region 
with the M(7r + 7r~) requirement relaxed; the histogram indicates the side-band determined 
background, which is parameterized by the fourth-order polynomial shown in the figure 
as a dot-dashed curve. The solid curve in Fig. El shows the result of a fit to the M(7r + 7r - ) 
distribution that uses the background function plus an acceptance- weighted p BW line-shape 
with an S- wave cut-off factor at the upper kinematic boundary 1J ; the dashed curve shows 
the fit with a P-wave cut-off factor. The .S'-wave case fits the data well: x 2, /d.o.f. = 43.1/39 
(CL=28%). The P-wave fit is much poorer, ^jd.o.j. = 71.0/39 (CL=0.1%), indicating 
that J ++ is strongly favored over J~ + . 

In summary, we find that with reasonable assumptions and a sample of 47 X — > 7i + 7i^J/ip 
signal events, we can rule out the J PC = 0~ + and ++ assignments for the X(3872) based on 
angular correlations among the final state particles. In addition, the M(tt + 7t~) distribution 
is inconsistent with all J~ + assignments. 

The results reported here, taken together with the observation of the X(3872) — > 'jJ/i/j 
decay mode 0, rule out all J PC assignments with J < 2 other than 1 ++ and 2 ++ . The decay 
angular distributions and 7r + 7r~ invariant mass distribution agree well with expectations for 
the the 1 ++ assignment. The 2 ++ assignment is not seriously challenged by any of the tests 
reported here, but is made rather unlikely by Belle's recently reported evidence for the decay 
X(3872) — ► D°D 7i° |l5|. The formation of 2 ++ from three pseudoscalars requires at least 
one combination to be in a P-wave. Thus, the near-threshold production of D°D°7i° would 
be suppressed by an £ = 2 centrifugal barrier. 

The 1 ++ charmonium x'd state is an unlikely assignment for the X(3872). Potential 
model predictions for the x'& m ass range from 3953 MeV ~ 3990 MeV [16], well above the 
X(3872) mass. The potential model masses are expected to be modified by coupling to 
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open-charm states. A coupled-channel calculation of open-charm-induced splittings for the 
x'd yields an upward mass shift of +28 MeV jl7l ]. 

The decay x'd ~~ > n + n~J/ip would proceed via pJ/ijj and violate isospin. The only 
well established isospin-violating hadronic transition in the charmonium system is ip(2S)_—> 
TT°J/ip, which has a measured partial width of T(ip(2S) — > n°J/tp) = 0.27 ± 0.06 keV [l2|. 
This is small compared to the expected total width of an M = 3872 MeV x'& °f m ore than 
1 MeV 0Q. A decay mode with a partial width this small would thus have a branching 
fraction that is less than 0.1%. This contradicts the recent BaBar 90% confidence lower 
limit of £(AT(3872) -»■ n+n-J/ip) > 4.3% EL Godfrey and Barnes calculate a partial 
width for an M = 3872 MeV x'd to be 11 KeV jig , more than an order-of- magnitude larger 
than that for the isospin violating ip(2S) — > 7r°J/"0 transition. Thus, one expects the jJ/ip 
decay to be stronger than pJ/ip. This is contradicted by our measurement: r(X(3872) — > 
7J/V>)/r(V(3872) -> Ti+Ti- J/tp) = 0.14 ± 0.05 Q. 

The 1 ++ a ssig nment is favored by models that treat the X(3872) as a molecule-like D°D*° 
bound state |l9l. Eo|. These models predict strong isospin violations and a "fJ/ip branching 
fraction that is much less than that for n+n-J/ip |2JJ, in agreement with observations. 
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